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ABSTRACT 

We present results of survey of interstellar absorptions towards supernova remnant 
(SNR) RX J0852.0-4622. The distribution of Ki absorbers along the distance of the 
background stars is indicative of a local region (d < 600 pc) strongly depopulated 
by K I line-absorbing clouds. This fact is supported by the behavior of the interstellar 
extinction. We find four high- velocity Call components with velocities of > 100 km s _1 
towards three stars and identify them with shocked clouds of Vela SNR. We reveal and 
measure acceleration of two shocked clouds at the approaching and receding sides of 
Vela SNR along the same sight line. The clouds acceleration, velocity, and Ca II column 
density are used to probe cloud parameters. The total hydrogen column density of 
both accelerating clouds is found to be similar (~6x 10 17 cm~ 2 ) which indicates that 
possibly there is a significant amount of small-size clouds in the vicinity of Vela SNR. 
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1 INTRODUCTION 

Optical absorption spectroscopy is a classical tool to probe 
interstellar (IS) clouds, viz., velocities, Doppler broadening, 
column density o f atomic and molecular s pecies, and tem - 
poral variations l|Danks fc Sembachl Il995l ; ICrawford| [2003). 
The direction of Vela SNR is unique among other regions of 
sky in one important respect: apart from usual low- velocity 
absorption components we see th ere fast-moving absorbers 
with velocities 100-200 km s" 1 (IWarierstein fc SilkJ [l97ll ; 
iDanks fc Sembachlll995l ; ICha fc Sembachl2000r i. These com- 
ponents are identified with IS clouds shocked in the blast 
wave of Vela SNR. An exciting feature of these high- velocity 
components is their temporal variability (jPanks fc Sembachl 
11993 : ICha fc Sembachl [2000). Among these variations of 
particular interest is the evidence of cloud acceleration 
l|Cha fc Sembachll2000h which probably reflects ongoing pro- 
cess of blast wave/cloud interaction in Vela SNR. Surpris- 
ingly enough, no attempts have been made so far to utilise 
acceleration of shocked clouds as an important observational 
constraint of cloud parameters. 

Recently, we performed low- and high-resolution spec- 
troscopic survey of 14 stars in the field of RX J0852. 0-4622 
aka Vela Jr. The latter is superimposed on the south-east 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, 080.D-0012(A) 



part of the Vela SNR. The analysis of low-resolution ob- 
servations in the re gion 3700-4000 A was reported earlier 
l|lvudin et al.ll2010r i. Here we address the high resolution 
spectral data which cover much broader spectral range. We 
find the results of analysis of these spectra are of great in- 
terest by several reasons. Firstly, our survey covers poorly 
observed region of Vela SNR. Onl y two stars HP 7530 9 and 
HD 75821 of the previous survey (ICha fc Sembachll200Ch fall 
into this region. Secondly, our spectra cover essentially larger 
spectral range than any other previous survey in the Vela 
direction. This permits us to obtain more extensive picture 
of interstellar atomic and molecular absorbers in the Vela 
SNR direction. Thirdly, an d most importan t ly, H D 75309, 
the common star of our and lCha fc Sembachl (|2000h survey, 
shows high- velocity Call absorption at -120 km s" 1 and 
+ 120 km s _1 presumably associat ed with shocked c louds . 
Our data combined with those of ICha fc Sembachl (|2000h 
survey cover large time interval of 15 yr which open up a 
possibility to measure rather precisely the velocity change 
of shocked clouds, and thus to recover the clouds acceler- 
ation. In turn, the cloud acceleration value in combination 
with the velocity and Call column density provides us with 
a new diagnostic tool for probing cloud parameters. 

The paper is organized as follows. We first describe ob- 
servations and the extraction of interstellar lines (Section^. 
In two cases the stars turn out slow rotators and, therefore, 
detailed stellar atmosphere analysis is used to recover inter- 
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Table 1. List of observed stars 



Star RA Peel m v SpType (B -V) A v d S/N 

(eq. 2000.0) (mag) (mag) (mag) (pc) all Call Nal ET 



HD 75309 


08 47 28.0 


-46 27 04 


7.84 


B2Ib/II 


0.01 


0.8 


1900±300 


113 


38 


146 


114 


HD 75820 


08 50 26.0 


-46 14 53 


8.64 


B9V 


-0.02 


0.2 


470±100 


79 


30 


101 


78 


HD 75873 


08 50 48.8 


-46 18 36 


8.10 


A3II/III 


0.38 


1.2 


1400±200 


67 


18 


88 


82 


HD 75955 


08 51 26.0 


-45 37 23 


7.73 


B9V 


-0.01 


0.2 


320± 70 


65 


20 


83 


61 


HD 75968 


08 51 32.8 


-46 36 36 


8.14 


B9III/IV 


-0.12 


0.0 


570±140 


45 


14 


61 


57 


HD 76060 


08 52 02.4 


-46 17 20 


7.88 


B8IV/V 


-0.09 


0.1 


390± 90 


101 


34 


130 


97 


HD 76589 


08 55 23.0 


-46 53 28 


8.34 


B9IV 


-0.05 


0.1 


390± 90 


67 


17 


87 


66 


HD 76649 


08 55 50.4 


-46 20 30 


8.33 


B7II/III 


0.14 


0.8 


640±110 


57 


14 


74 


65 


HD 76744 


08 56 18.2 


-46 19 57 


8.69 


AOV 


0.08 


0.5 


270± 50 


61 


18 


78 


59 


CD-454590 


08 49 35.5 


-46 23 18 


9.58 


B5 


0.20 


1.3 


2400±300 


11 


11 


53 


13 


CD-454606 


08 50 15.0 


-45 31 22 


8.96 


B0.5V 


0.38 


2.0 


1670±160 


63 


18 


82 


72 


CD-454645 


08 51 34.9 


-46 09 54 


10.32 


AO 


0.20 


0.4 


330± 70 


38 


11 


49 


41 


CD-454676 


08 53 22.0 


-46 02 09 


8.93 


B0.5III 


0.77 


3.2 


1080±150 


38 


13 


52 


54 


CD-464666 


08 50 44.3 


-46 38 11 


9.81 


AMI 


0.60 


2.1 


5700±500 


37 


12 


50 


51 



stellar absorption. The parameters of interstellar absorption 
components and their analysis are presented in Section [3] In 
Section [4] we address the shocked clouds in the sight line of 
HD 75309. The cloud acceleration, velocity and Call column 
density are used for probing cloud parameters on the bases 
of cloud shock model in which we include all the relevant 
physics. 



2 OBSERVATIONS AND DATA REDUCTION 

The spect ra of 14 stars in di rection of Vela Jr. [for star po- 
sition see i|lvudin et al.ll2010h ] were obtained in March 2008 
on the 3.6-m ESO NTT telescope equipped with EMMI 
spectrograph (echelle grating #14, cross-dispersion grat- 
ing #3, CCD mosaic 2076x4110). The full coverage is 3850- 
8620 A with the gap at 5570-5650 A; the central dispersion 
is 0.02 A/pix, the resolving power is R = A/AA=88 600 
(Au=3.4 km s- 1 ). In Table Q] stars are listed with their 
HD/CD name, equatorial coordinates, ^/-magnitude, spec- 
tral type, B — V color, and signal-to-noise ratio (S/N), av- 
erage and at Can 3933,3968 A, Nai 5889,5895 A, and Ki 
7664,7698 A positions. Preliminary processing of CCD im- 
ages, spectra extracting, wavelength calibration with ThAr- 
lamp, normalizing and merging of orders is performed using 
echelle package of MIDAS. In order to normalize the spec- 
tra a blaze function of the echelle is recovered using a cali- 
brated flux of the standard star HD 6075 3. Its stellar fl ux is 
synthesiz ed by applying A TLA S9 code (|Kuruczlll993h and 
SynthVb (|Tsvmbal et al.ll2003l ). A stellar atmosphere model 
of standard star is calculated for the solar metallicity with 
parameters T e // = 16200K, log q=3.56, V^sin ^=24 km s~ , 
and interstellar absorption Av = 0.25 jlvudin et al.ll2010l ). 

For most stars we used stellar synthetic spectrum as- 
suming solar element abundances because these objects be- 
long to the thin disk of Galaxy. We find that small devi- 
ations (< 0.3 dex) from the solar abundance do not affect 
the extracted interstellar lines, if the stellar rotation is high 
(> 30 km s" 1 ). In case of star HD 75873 the rotation is slow, 
Vsin i— 3 km s _1 , and, therefore, a stellar synthetic spec- 
trum with carefully determined chemical composition should 
be used to extract interstellar lines. A similar procedure is 
applied to CD-464666. 



The stellar parameters of HD 75873 are found by em- 
ploying a standard requirement that the chemical compo- 
sition should be stable against variations of the line exci- 
tation potential and ionization state. The ATLAS9 code 
is used to derive parame ters of the stellar atmosphere and 
WIDTH9 jKuruczll 19931 ) for the abundance determination. 
The stellar parameters infered in this way (T e //=8900K, 
log g=2.50) are consistent with those determi ned earlier us- 
ing only hydrogen lines jlvudin et al.l |2010| ). The micro- 
turbulent velocity (Vt=2.5 km s" 1 ) is derived using re- 
quirement that the Fe/H ratio should not depend on the 
intensity of Fell lines. The recovered element abundance 
slightly exceeds the solar value; the Fe overabundance is 
+0.08±0.11 dex. In case of CD-464666, most distant star 
of the program, we find T e// = 10100±300 K, logg=2.0±0.2, 
and Vt=2.5±0.5 km s~ . The distance we infer using the 
spectral parallax method is 5800=1=9 00 pc, in agreement with 
the previous estimate 5700±500 pc (|lvudin et alj|2010l ). The 
inferred element abundance is close to the solar one; Fe over- 
abundance is of +0.05±0.11 dex. 

The extracted spectrum of the interstellar absorption 
(Fig. [TJ is represented by a residual flux 

F{\) = 1 - (F syn {\) - F obs (A))/F cont (A) , 

where Fobs is the observed normalized flux of the stellar 
spectrum, F syn is the model normalized flux calculated usin g 
SynthV package and VALD line list (jKupka et al.lll999r i. 
F cont is the continuum flux. We extracted telluric lines in 
case of Nal 5895 A, K I 7664 A, and diffuse interstellar bands 
(DIBs). 



3 DATA ANALYSIS 

3.1 Overview of IS absorptions 

The interstellar spectra show absorption lines of Call 3933, 
3968 A, Nai 5889, 5895 A, Ki 7664, 7698 A, CH 4300 A, 
CH+ 4232 A, and DIBs. The Cal 4226 A line is not de- 
tected, in accord w ith the well known weakness of this line 
ijWeltv et al.ll2003l ). In Fig. [2] we present a zoom of Nal, 
Can, Ki, CH, and CH + interstellar lines for the interest- 
ing case of HD 75309 (d = 1900 pc) showing high- velocity 
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Figure 1. Interstellar absorption lines of Kl, Call, and Nal in spectra of program stars. The star distance increases from left to right 
and from top to bottom. In the right-bottom corner we show position of spectral lines. All continuum fluxes are normalized to unity. 
Spectra are shifted in vertical direction for convenience. 
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Table 2. Equivalent widths (mA) of interstellar atomic and molecular lines 



otar 


Can 


Nal 


IV 1 




Uxl 




3933 A 


3968 A 


5889 A 


5895 A 


(004 A 


fbyo A 


a Qnn A 
4oUU A 




HD 75309 


297 


131 


567 


450 


57 


52 


6 


7 


HD 75820 


96 


16 


110 


50 








_ 


_ 


HD 75873 


101 


100 


536 


454 


204 


161 


26 


10 


HD 75955 


87 


21 


199 


116 










HD 75968 


62 


17 


122 


74 










HD 76060 


108 


22 


192 


110 










HD 76589 


85 


28 


115 


63 










HD 76649 


41 


65 


383 


283 


82 


89 


15 


8 


HD 76744 


77 




101 


50 










CD-454590 


431 


113 


568 


533 


164 


132 


22 


20 


CD-454606 


368 


189 


632 


496 


198 


153 


22 


11 


CD-454645 


17 


10 


244 


169 










CD-454676 


156 


107 


470 


427 


228 


208 


37 


6 


CD-464666 


283 


147 


780 


694 


177 


136 


6 


17 



Table 3. Equivalent widths (mA) of DIBs 



Star 


5780 A 


5797 A 


5850 A 


6196 A 


6203 A 


6270 A 


6284 A 


6376 A 


6379 A 


6614 A 


6660 A 


7224 A 


HD 75309 


143 


20 




11 


43 


26 


308 




16 


52 


10 


31 


HD 75820 


























HD 75873 


321 


75 


25 


29 


73 


62 


598 


16 


41 


118 


24 


56 


HD 75955 




12 






















HD 75968 


23 












27 












HD 76060 


22 












24 












HD 76589 


18 
























HD 76649 


114 


27 






19 




253 




11 








HD 76744 


10 
























CD-454590 


246 


97 


IS 


27 


62 


47 


620 


15 


54 


107 


30 


65 


CD-454606 


261 


71 


21 


24 


57 


41 


409 


18 


46 


128 


21 


14 


CD-454645 














63 












CD-454676 


334 


131 


35 


30 


82 


37 


769 


20 


71 


151 


21 


86 


CD-464666 


485 


125 


30 


35 


113 


59 


966 


23 


58 


161 


33 


94 



(\v\ > 100 km s _1 ) Call components that will be addressed 
below. Total equivalent widths of the detected interstellar 
absorptions are given in Table [2] 

The interstellar Nal doublet is detected in the spectra 
of all stars (Fig.[T}. This line, as expected, becomes stronger 
as star distance increases. In most cases Call doublet is sig- 
nificantly weaker than that of Na I which is well known fact 
and usually explained as an effect of a strong Ca deple- 
tion onto dust grains. Kl doublet is seen only in stars with 
large extinction and strong molecular absorptions of CH and 
CH+ . The equivalent width of K I 7698 A absorption tightly 
correlates with the total equivalent width of molecular lines 
of CH and CH + . (Table [2}. This correlation agrees with 
the strong correlation between Kl and CN column densities 
l|Weltv fc Hobbsll200ll ). 

The spectra arranged according to the star distance 
in Fig. [T] reveal a striking behavior of Kl absorptions: 
they are not detected for stars with distances d < 600 pc 
(iV(Ki) < 5 x 10 10 cm -2 ) and then sharply get strong for 
stars with distances d > 600 pc (Fig. [3]a) . The "jump" in 
the distribution of K I absorbers along the distance is indica- 
tive of a "cavity" in the region d < 600 pc underpopulated 
by cool clouds. This conjecture is supported by a similar- 



behavior of the interstellar extinction although in this case 
the threshold effect is not as strong as for Kl absorbers 
(Fig. Ob). 

Twelve diffuse interstellar bands are detected in our 
sample and listed in Table [3] Empty place in the ta- 
ble indicates that this particular DIB is absent or its 
equivalent width is smaller than 10 mA. Most DIBs with 
A <5500 A cannot be extracted because of low S/N ratio. 
In an extensive surve y of DIBs in spectra of ~ 200 stars 
|Friedman et al.ll201ll ) finds a strong correlation of equiva- 
lent widths (W in mA) of 6284 A and 5780 A DIBs (cor- 
relation coefficient r = 0.96) with the linear fit W^(6284) = 
a + (5780), where a = 28. 24 ±5.8 and b = 2.32 ±0.03. In 
our small sample the correlation even stronger (r = 0.99), 
while the parameters of the linear fit are a bit different: 
a = 1.5 ± 46 and b = 2.03 ± 0.18. Yet the relatively small 
value of a term in our case seems to be more easy to con- 
form with the expected zero value of W in case of vanishing 
molecular abundance. 

3.2 Component parameters 

ICrawfordl (|200lh demonstrated upon the bases of high- 
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Figure 2. Interstellar absorption lines in spectra of HD 75309. 
All continuum fluxes are normalized to unity. Spectra are shifted 
in vertical direction for convenience. 
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Figure 4. Velocity distribution of resolved components of inter- 
stellar absorptions of Kl, Nal, and Call. 
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Figure 3. Equivalent width of Kl absorption (upper panel) and 
extinction (lower panel) vs. star distance. Arrows in upper panel 
show upper limits. 

resolution (Ad «0.3 km s _1 ) observations that Doppler pa- 
rameters of interstellar line components lie in the range of 
0.3 ^ b ^ 2.2 km s _1 . Our spectra have lower resolution 
(Av — 3.4 km s _1 ) so most of the interstellar line com- 
ponents are expected to be unresolved. This precludes re- 
liable determination of the Doppler parameter of a single 
line. Fortunately, in cases of Nal, Can, and Kl we have in 
hand a doublet ratio; this provides us with the additional 
relation which permits one to recover both TV and 6 even for 
unresolved lines. Generally, the blend of M components is 
decomposed using standard relations 

F(v) = (e~ r W * 9) (1) 

M / _ \ 

t(v) = ff0 /A£jV^(^-^J , (2) 

where the first equation is a convolution of the line pro- 
file with the instrumental profile O. In the second equa- 
tion <7o = 0.0265 is the integrated cross-section parame- 



ter, / is the oscilla tor strength taken from VA LD database 
|Kupka et all 1 19991 ) for atomic lines and from lGredel et all 
jl993h for molecular lines, Nk is the column density of k- 
th component, and <j>{u/b) = (1/b^/w) exp(— u 2 /6 2 ) is the 
Gaussian velocity distribution of an absorbing gas in a cloud 
along the sight line. 

Parameters of interstellar absorption components de- 
rived by the minimization procedure are listed in Table [4] 
The radial velocity is reduced to the local standard of rest 
(LSR) using t he Sun velocity to wards Vela Jr. Avi SI — 
-15.7 km s -1 dFamaev et al.ll2005h . For each atomic species 
three columns (v\ BT , b, and N) are given. The Doppler pa- 
rameter for molecular lines is omitted because it cannot be 
reliably determined. Data with the large errors are marked 
by the colon. The typical vi SI error is 0.1-0.2 km s _1 (0.5- 
1 km s" 1 in case of the low accuracy), b error is 0.1- 
0.4 km s~ (1-1.5 km s _1 in case of the low accuracy); for 
N a typical relative error is 10% and 30% in case of the 
low accuracy. The blends with unresolved saturated lines 
are marked by italic with the estimated velocities shown in 
parentheses. 

The average Nal column density of components in 
our sample (3.5 x 10 11 cm" 2 ) is larger compared to 1.2 x 
10 11 c m" 2 , the ave r age v alue for the high-resolution sur- 
vey of IWeltv et all ()l994h . We attribute this disparity to 
the low resolution of our spectra which leads to the omis- 
sion of the significant number of thin cloud components in 
absorption blends. This effect is apparent also in the ve- 
locity distribution (Fig. [4} . K I absorptions are rather weak 
so the components are recovered confidently. The velocity 
dispersion of Kl absorbers a v w 8 km s" 1 is the same as 
the value 7.8 km s" 1 found for an e xtensive survey for K l 
components with logA r < 11.3 cm" 2 (jWeltv fc HobbJ200ll ). 
For stronger components these authors find somewhat lower 
value a v = 5.8 km s~ . Our small sample precludes more re- 
fined analysis of the velocity distribution of strong K I com- 
ponents. In case of Na I the velocity distribution (Fig. [4| is 
rather broad with a v ~ 18 km s~\ twic e as la rge compared 
to a v ss 9 km s _1 found bv lWeltv et all [| 19941 ) in an exten- 
sive survey of the Nal IS absorption. The striking disparity 
is caused probably by the relatively low resolution and S /N 
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Table 4. Parameters of the resolved components of interstellar absorptions 



Star 


Can 






Nal 






Ki 




CH 




CH+ 






v lsI b 


Nil 




b 


N11 v lsr 


b 


N11 




N11 


flsr 


N11 




(kms^ 1 ) (km a -1 


) (cm" 2 ) 


(km s- 1 ) 


(km a" 1 ) 


(cm 


2 ) (km a" 1 ) 


(km s- 1 ) 


(cm" 2 ) 


(km s- 1 ) 


cm" 2 ) 


(kms- 1 ) ( 


cm- 2 ) 


ilU (oouy 


-1Z ( . 1 1.0 


AX) 


-17: 






2.7 


1 ^ 
l.O 


1.0 


o.u. 


/ 0. 


D.D 






-77: 2: 


1: 


-2.6 






8.6 


2.0 


0.6 












-65: 1.5: 


0.5: 


2... 11 (5.0: 8.8:) 




















-28...25(-22: -9: 


2.5: 12.5:) 


20: 


3: 


3: 


















33... 55 


























123.9 1.5: 


0.6 






















HD 75820 


-14.9 


3.0: 


-27.1: 


1.0: 


0.5 


















2.0 


3.0 


-22.2: 


1.0: 


0.5 






















-13.6: 


1.5 


0.6 






















-7.6: 


1.0: 


0.9 






















-1.0 


1.1 


1.1 






















4.5 


1.0: 


0.5 






















8.2 


1.0: 


0.4 
















HD 75873 


1.4: 


13: 


-25.6 


3.0: 


0.6 


3.2 


1.5 


4.4 


4.5: 


400: 


6.1 


158 








-6.8 


2.3: 


4.8 


9.1 


1.5 


7.6 
















-S...25 






15.4 


0.7: 


1.1 
















30.2 


3.0: 


1.5: 














HD 75955 


-21: 


3.5 


-39.2 


1.4 


0.4: 
















-3: 


5.5 


-19.3 


1.3 


0.5: 




















-1.1 


1.4 


16.5 




















7.5 


1.7 


1.2 
















HD 75968 


1.0: 


6.0: 


-12.8 


1.0: 


0.4 






















-6.7 


1.5: 


0.4 






















-1.6 


1.0 


1.7 






















10.8 


1.5: 


0.7 






















27.2 


0.5: 


0.8 
















HD 76060 


-121.8 4: 


4: 


-0.8 


1.5 


5.6 


















-15...7(-8: -1:) 




9.1 


2.5: 


1.0: 
















26.4: 7: 


6.3: 


33.2 


0.9 


3.9 
















HD 76589 


-27.. .13 


-21.1 


2.0 


0.5 




















35: 


1.4: 


-7.5 


1.5 


0.5: 




















-1.2 


1.4 


2.6 






















3.8 


0.7 


0.8 






















16.1 


3.5 


0.7 
















HD 76649 






-2.2 


1.3 


3.8 


6.9 


1.0: 


6.7: 


8.1 


220 


8.0 


130 








1...16 






10.6 


0.3: 


1.5: 










HD 76744 


-5.0: 




-22. ..-15 


























-6.5 


1.9 


0.6 






















-1.7 


1.5 


1.2 






















5.2 


2.6 
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of our spectra which precludes decomposition of saturated 
blends of Na I absorptions. As a result we strongly underes- 
timate the number of low- velocity components of Nal. 

The large velocity dispersion of Can absorption com- 
ponents (Fig. [4j stems from two effects. The first is the 
same selection effect as in case of Nal absorption, while 
the second effect i s the well known phenomenon discussed 
iRoutlv fc Spitzed l|l952h : the larger ratio iV(Caii)/iV(Nai) 
for components with higher velocity. A standard explana- 
tion of the Routly-Spitzer effect is the dust destruction in 
the dynamically disturbed clouds which results in release of 
gaseous Ca. The fact that we see high-velocity components 
in Ca II lines and not in Na I lines indicates that high- velocity 
clouds have low neutral hydrogen column density in which 
case the equivalent width of Na I absorption turns out below 
the detection limit of our data. 

In some cases positive velocities of ui sr < 50 km s~ x are 
related with a velocity gradient caused by the galactic ro- 
tation. For the most distant star CD-464666 (d — 5700 pc) 
adopting galactocentric distance of 8.5 kpc and Galaxy ro- 
tational velocity of 220 km s _1 we find that Galactic rota- 
tion can be responsible for the positive LSR velocities in the 
range of v ^ 50 ±4 km s . Remarkably, this value is consis- 
tent with the upper limit of the observed velocity interval of 
the Can absorption (-25. ..+50 km s _1 ) towards CD-464666. 

Two nearest stars of our sample, CD-454645 (d — 
330 ± 70 pc) and HD 76060 (d = 335 ± 63 pc), show a 
strong Nal component with similar velocity (32.6 km s -1 
and 33.2 km s _1 ), Doppler parameter (1.1 km s _1 and 
0.9 km s -1 ), and column density (3.9 x 10 11 cm -2 ). This 
suggests that we see both stars through the same IS cloud. 
For the distance of ~330 pc the angular separation 9' be- 
tween these stars corresponds to the linear scale of 0.85 pc. 
We conclude therefore that the tangential size of this cloud 
> 1 pc. The large LSR velocity (33 km s _1 ) indicates that 
this cloud belongs to the dynamically disturbed ISM. Re- 
markably, however, that we do not detect similar component 
in Call suggested by Routly-Spitzer effect. This indicated 
that this fast-moving material contains significant amount 
dust with the most Ca locked in grains. 

Of particular interest are four fast-moving Ca II compo- 
nents (Table[3]) with velocities exceeding 100 km s in sight 
lines of three stars, specifically, -177 km s _1 (CD-454676), 
-125 and +120 km s _1 (HD 75309), and -122 km s _1 
(HD 76060). Following the earlie r conjecture invoked for 
similar high-velocity components JWallerstein fc Silklll97ll ; 
IJenkins et all Il976r , ICha fc Sembachl |2000| ) towards Vela 
SNR we attribute the fast-moving material to the IS clouds 
shocked by Vela SNR. 

Alternatively, high velocity clouds of our survey might 
be attributed to the SNR RX J0852. 0-4622. We consider, 
however, this possibility unlikely because RX J0852. 0-4622 
is a young SNR that expands with the veloci ty of ~ 
10 4 km s" 1 jlvudin et al ] |l998l ; lAschenbachlll998l ). In this 
case velocities of shocked clouds unlikely would be the same 
as in old Vela SNR. The identification of detected high- 
velocity components with Vela SNR is thus preferred. 
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Table 5. Evolution of high-velocity components in spectra of 
HD 75309 



Year 




cloud A 




cloud B 






W 


N 


«->lsr 


W 


N 




(km s -1 ) 


(mA) 


(lO^cm" 2 ) 


(km s" 1 ) 


(mA) 


(10 n cm- 2 ) 


1993 


-121.9+0.3 


12+1 


1.4+0.1 


119.3±0.3 


9+1 


1.0+0.1 


1996 


-123.8±0.1 


16±1 


1.8+0.1 


119.9+0.3 


7+1 


0.8±0.1 


2008 


-127.1+0.1 


16±1 


1.8+0.1 


123.9+0.1 


5+1 


0.6±0.1 



T 




-140 -120 -100 100 120 140 

V, , kms V, , kms 

1st lsr 

Figure 5. The profiles of Can absorption in spectra of star 
HD 75309 for three epochs. For all the epochs Call K are shown. 
In the inset we show both K (thick) and H Call lines in 2008 
spectrum to demonstrate that the noise does not affect the nar- 
row component position. 



3.3 Acceleration of high-velocity clouds in sight 
line of HD 75309 

Interstellar lines towards HD 75309 were o bserved in Febru- 
ary 1993 and January/February 1996 (|Cha fc Sembachl 
2000) with the spectral resolution comparable with that of 
our data. Authors noticed increase of a column density of 
Call fast component with velocity -120 km s~ but did not 
mention any velocity variation. We combine their data with 
our observations to look at the velocity change on the time 
sc ale of 12-15 yr. It s hould be noted that velocities reported 
by ICha fc Sembachl (|2000h are larger compared to ours by 
4 km s _1 for all the components. The major part of this 
difference (2.6 km s -1 ) is explained by the new value we use 
for the solar motion l|Famaev et al.|[2005h . The remaining 

1.4 km s _1 which is comparable with one pixel is unclear but 
could be related with the wavelength calibration uncertainty 
of both spectra. Despite this uncertainty the velocity change 
of the fast components can be measured with high precision 
via the cross-correlation technique because slow components 
of interstellar absorption lines remain unchanged. This pro- 
cedure permits us to establish a common zero point in the 
velocity space and thus to measure the velocity change over 
the past period of time. 

The high-velocity components of -127 km s" 1 and 
+ 124 km s" 1 a ccordi ng to our data and to those of 
ICha fc Semblchl l|2000l) with the common zero point are 
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shown in Fig [5] We name the approaching component 
"cloud A" and receding component "cloud B". The plots 
clearly demonstrate that both components accelerate. Re- 
markably, in the case of cloud A both Can K and H lines 
in the 2008 spectrum show similar narrow component thus 
strengthening the reliability of the velocity determination. 
Using Gaussian fitting we measure velocity of both clouds 
at three epochs (Table [5]). In the table we give also W and 
A; the dat a on W and A f or ep ochs 1993 and 1995 are 
taken from ICha fc Sembacbl ([2000). The average accelera- 
tion between 1993 and 2008 is -0.33±0.06 km s" 1 yr" 1 for 
cloud A and +0.38±0.06 km s^yr" 1 for cloud B. The ap- 
proaching component seems to change its width while the 
equivalent width is preserved. The equivalent width of the 
receding component has noticeably decreased. 

The impact parameter (in SNR radius) of the sight line 
of HD 75309 relative to the Vela SNR center is p/R ~ 0.5. 
Adopting p/R = 0.5 and assuming the model of a plane slab 
perpendicular to the SNR radius we find for the cloud A the 
deprojected velocity v c = 147 km s _1 , Can column den- 
sity .A(Can) = 1.6 x 10 11 cm" 2 , and the cloud accelera- 
tion g=0.38 km s" 1 yr" 1 (or 1.2 x 10~ 3 cm s~ 2 ). For the 
cloud B these values are v c = 142 km s~ , A(Can) = 6.9 x 
10 10 cm -2 , and g = 0.44 km s _1 yr -1 (or 1.4x 10~ 3 cm s -2 ). 
The deprojected values will be used below for the constrain- 
ing the cloud parameters. 



4 PARAMETERS OF SHOCKED CLOUDS 

The question we address here is what are the cloud param- 
eters (density and size) which might explain simultaneously 
the Can column density and the acceleration of shocked 
clouds. Below we rely on the assumption that the interstel- 
lar clouds are shocked and accelerated by the blast wave of 
Vela SNR. We start with the general picture and then use a 
model of the cloud shock to constrain the cloud parameters. 



4.1 General picture 

We follow a view that the Vela SNR is the result of an ex- 
pansion of the Sedov blast wave in the cloudy IS medium 
jBocchino et all 1 19991 ). In our estimates and modelling we 
approximate sub-parsec clouds by a homogeneous sphere 
which is the idealization for the clumpy cloud. The assump- 
tion of a spherical cloud is usually adopted for the an alysis 
of a blast wave/cloud interaction (|Klein et al.l lT994). par- 
ticularly in case of a sh ocked cloud in the Cygnus Loop 
llPatnaude fc Fesenll2005h . In fact, since we are interested 
primarily in the cloud size along the blast wave propaga- 
tion we do not concern about the clump shape and tangen- 
tial size. Some justification for the choice of clumpy shape 
provide 3- dimension al magneto-hydrodynamic simulations 
jKritsuk et all l201lh which demonstrate that initially ho- 
mogeneous medium evolves towards a two-phase structure in 
which cool phase grows filamentary but later on the clumpy 
structure gets dominant. Yet we admit that in some cases a 
cloud shape may be close to that of a filament or sheet. 

A blast wave propagating with the speed in an in- 
tercloud medium with the density pio dr ives a cloud shock 
wave with the speed v B « Vh (pio/pco) 1/2 l|Klein et al.lll99l ) 



where p c o is the undisturbed cloud density. (Note, we dis- 
tinguish between observed cloud velocity v c and the cloud 
shock veloci ty v s ). Adopting the blast wave speed = 
600 km s _1 |Bocchino et al.lll999l ) and the cloud shock ve- 
locity v B = 140 km s _1 one finds from the momentum 
conservation that the cloud-to-intercloud density contrast 
X = Pco/pio ~ ( vb/vs) 2 ~ 18. For a t ypical intercloud den- 
sity TliO = 0.1 jBocchino et al l Il999h the cloud density is 
then n c o ~ 2 cm -3 . 

The interaction of the blast wave with the cloud of 
the radius a can be divided roughly into three phases 
(Klcir Tet al.lll994l ). At the first stage the cloud shock prop- 
agates through the cloud; this stage ends up with the cloud 
crush in a "cloud crushing time" t cc — a/v s . At the sec- 
ond stage the shocked cloud is accelerated by the dynami- 
cal pressure of the blast wave. This is accompanied by the 
growth of the Kelvin-Helmholtz instability that leads to the 
cloud fragmentation. The third stage is the final stage of 
fragmentation which ends up with the full destruction of a 
cloud and mixing of cloud fragments with the blast wave 
flow. For the de nsity contrast 10 ^ \ ^ 100 the cloud life 
time is id ~ 4t cc l|Klein et al.| [l994). It should be noted that 
the shocked cloud may not show acceleration. This would be 
the case if the time t cc is so large that the cloud has enough 
time to cross the blast wave before it gets fragmented. The 
fact that we see the accelerating Ca II line-absorbing gas in- 
dicates that the cloud has been already completely shocked 
and still resides inside the dense part of the blast wave. We 
however are not able to exactly specify the stage at which 
the shocked cloud is detected. The conservative statement 
would be that the accelerating cloud is caught at the stage 
t cc < t < 4t cc after the interaction between the blast wave 
and the cloud has turned on. 

The cloud acceleration combined with an additional in- 
formation about Vela SNR can provide us with an estimate 
of the cloud column density via the equation of a cloud mo- 
tion driven by the blast wave. For the blast wave postshock 
velocity of (3/4)t>b and the postshock density of 4riio the 
equation of the cloud motion reads 

gN H = 4n i0 f ~«b - «cl , (3) 

where Ah is the cloud total hydrogen column density 
normal to the shock plane. Inserting into this equation 
g » 10~ 3 cm s" 2 , Vc = 140 km s" 1 , n i0 = 0.1 cm -3 , 
Vb = 600 km s -1 one infers the cloud column density 
Ah ~ 5 x 10 17 cm -2 . This exciting result shows that the 
accelerated cloud is unlike a typical Nal absorber which is 
associated with a c loud column density Ah ~ 8 x 10 19 cm -2 
jWeltv et al.lll994l ). The derived estimate of Ah combined 
with the cloud density n c o ~ 2 cm -3 implies the cloud ra- 
dius a = (3/4)AH/n c o ~ 2 x 10 17 cm assuming initial spher- 
ical cloud shape. For this radius the cloud crushing time is 
tcc — a/v c ~ 500 yr. 

The above estimates presume a constant dynamical 
pressure pdyn = pv 2 of the blast wave. This assumption, 
however, could easily break down because the dynamical 
pressure in the b last wave drops dow nstream rather steeply, 
p dyn oc (r/R) 11 l|McKee et al.lll987f ). Already at r = 0.977? 
the pressure is « 30% lower compared to the maximum 
value at the shock front. It is reasonable to adopt the post- 
shock layer A r ~ 0.03i? to be the width of the blast wave 
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in which the dynamical pressure is approximately constant. 
Because v c <C Vb, it takes At ~ Ar/vb ~800 yr for the 
cloud in the blast wave frame to cross the blast wave width 
A r = 0.037? « 0.5 pc, where we adopt R = 16 pc (C ha et all 
1 19991 ). The found time scale is only slightly larger than the 
estimated cloud crushing time t cc ~ 500 yr. This has an im- 
portant implication for the cloud sizes: the initial radius of 
accelerating clouds cannot significantly exceed 2 x 10 17 cm; 
otherwise the cloud crosses the blast wave before it gets 
shocked. In this regard it might well be that not all the 
high-velocity shocked clouds accelerate; some of them might 
be large enough to cross the blast wave before they get frag- 
mented. These clouds obviously should not experience no- 
ticeable acceleration. 

4.2 Model outline 

The cloud postshock flow will be modelled in a steady 
plane shock approximation. The major parameters of the 
shock are the observed cloud shock velocity (v s ), preshock 
cloud density (n c o), preshock magnetic field (So), and 
preshock hydrogen ionization fraction (xq). The postshock 
flow then is determined by a standard set of equatio ns of 
mass, momentum and energy conservation (Cox 1972). We 
ad opt the parallel frozen- i n mag netic field. Cooling function 
of ISutherland fc Dopital (119931 ^ a nd the low-temperature 
(logT < 3.9 K) cooling function of lJura fc Dalgarnol (|l972T ) 
are used for the cooling rate calculation. We ignore effects 
of elements depletion onto dust in the cooling function be- 
cause in the temperature range of interest the cooling func- 
tion is determined by elements (H, He, C, N, O) which are 
not affected by the depletion. The hydrogen, helium, Ca, 
and Na ionization fractions are calculated by solving time- 
dependent kinetic equations which include radiative recom- 
bination, ionization by electron collisions, ionization by the 
interstellar ultraviolet radiation, by cosmic rays and back- 
ground X-rays. Among the latter two mechanisms the ion- 
ization by background X-rays dominates for small size clouds 
10 18 cm). We adopt the ioniz ation rate by backg round 
X-rays to be ( « 7 x 10" 16 s" 1 (jWolfire et al.lll995l ). The 
ionization of Ca n by the recombination Lyc radiation is also 
taken into account although this process is of minor impor- 
tance. Four Ca ions and three Na ions are included in the 
ionization calculations of these elements. 

The interstellar calcium is usually strongly de- 
pleted onto dust. We adopt depletion factors 5 = 
(Ca /H)/(Ca/H) ? = 2 x 1C T 4 for Ca and 5 = 0.1 for 
Na (I Savage fc Sembachlll996l ). The dust destruction in the 
shock is treated following standard recipes (|Tielens et al.l 
1 19941 ). We take into account the grain acceleration by the 
betatron mechanism (jshufi|[l977]) and the grain deceleration 
due to the gas-grain collision. The efficiency of the betatron 
acceleration is set to be unity, if the grain Larmor radius 
is small, r_e(dln B/dz) < 1, and zero otherwise. The ma- 
jor mechanisms of a dust destruction are the grain evap- 
oration in grain-grain collisions, nonthermal grain sputter- 
ing by collisions with He, and the thermal d ust sputtering 
jDraine fc Salpetedll979l ; iTielens et al.lll994l ). The adopted 
critical velocity of the vap orization in the grain-gra in col- 
lision is d„ = 30 km s" 1 (|Borkowski fc Dwekl fl995h . The 
surface binding energy relevant to the grain sputtering is 
set to be 4 eV, the value typical for amorphous carbon and 
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Figure 6. The postshock flow for the shock speed of 140 km s _1 
as a function of the downstream column density. Shown are (a) 
temperature (thick line) and density; (b) Ca depletion factor; 
(c) Can column density (thick solid), Nal column density (thin 
solid), observed Call column density (dotted), and Nal column 
density 3<r upper limit (triangle); (d) ionization fraction of H 
(thick solid), Call (thin solid), and Cal (dotted). 
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Figure 7. The same as Fig. \E\ but for the shock speed of 
100 km s" 1 

iron grains l)Tielens et al.lll994l ). Calcium locked in the dust 
is assumed to be homogeneously mixed with the grain ma- 
terial. 

An approximation of a single grain size with the initial 
grain radius of 10~ 5 cm is used. The preshock hydrogen 
ionization fraction is set to xo = 0.1 and the magnetic field 
is set to Bo = 3 x 10 -6 G; in fact, results are not very 
much sensitive to these parameters. In order to estimate the 
cloud acceleration we adopt the dynamical pressure of the 
blast wave to be p c v^ and the intercloud blast wave speed 
of 600 km s" 1 . 

4.3 Results of modelling 

We focus on the approaching cloud (cloud A) and consider 
two cases (I) the cloud was shocked recently with the shock 
speed of 140 km s , (II) the cloud was shocked a long time 
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ago by the cloud shock with the speed of 100 km s~ and 
afterwards was accelerated up to 147 km s _1 . We have ex- 
plored an extended parameter space and found an optimal 
models which meet constraints from both the Call column 
density and acceleration for the case I (Fig. [6]). and the 
case II (Fig. 0. 

In case I the observed velocity 147 km s -1 is reached 
in w 20 yr after the cloud crushing. The optimal value for 
the cloud preshock density is n c o = 4 cm -3 , while the in- 
tercloud density turns out to be nio = 0.22 cm -3 assuming 
Vb = 600 km s -1 . Shown in Fig. \6\ as functions of the to- 
tal hydrogen column density are the temperature, the total 
hydrogen number density, the Ca depletion parameter, the 
gaseous Call and Nal downstream column density, and the 
ionization fractions of H, Ca II, and Ca I. The column density 
of Call is sampled primarily in the cooling zone, at about 
JVh ~ 4 x 10 17 cm -2 . The rapid increase of the Call density 
at this distance is an outcome of the density increase in the 
cooling zone which is followed by the recombination of Ca ill 
into Call and the grain destruction. The dominant mecha- 
nism of the dust destruction is the non-thermal sputtering 
due to grain collisions with He atoms and ions. We find that 
results are robust with respect to the initial depletion fac- 
tor. The observed Call column density for the cloud A is at- 
tained in case I at Nn ~ 6.6 x 10 17 cm' 2 (Fig.[6jj). This can 
be considered as an estimate of the initial column density 
of the undisturbed cloud A. The acceleration of the shocked 
cloud with this column density is g = 1.2 x 10 -3 cm s -2 
which coincides within errors with the observed deprojected 
acceleration of the cloud A. For the adopted cloud column 
density the radius of the unshocked cloud is a = 1.2x 10 17 cm 
and the cloud crushing time is tcc = 370 yr. Remarkably, 
the model Na I column density is consistent with the 3<r up- 
per limit derived from the non-detection of the high- velocity 
Nal D2 absorption (Fig. EJ;). 

In the case II (Fig. [7} the observed Ca 11 column density 
and acceleration are satisfied for the total hydrogen column 
density of Nn = 5.5 x 10 17 cm -2 ; this value is presumably 
the column density of cloud A in case II. The cloud number 
density in this model is n c o = 6 cm -3 and the intercloud 
density is nio = 0.16 cm -3 assuming Vb = 600 km s _1 . 
The cloud is crushed at t cc = 290 yr while the observed 
velocity v c = 147 km s -1 is attained at t = lAt cc ~ 400 yr. 
Remarkably, both cases I and II suggest that the shocked 
cloud is observed at about 400 yr after the interaction began. 

The above models demonstrate that for the shock ve- 
locity in the range of 100-140 km s _1 the cloud total hy- 
drogen column density turns out to be in a narrow range 
of (5.5-6.6) x 10 17 cm -2 . The number density in the undis- 
turbed cloud is in the range of 4-6 cm~ 3 with the lower 
value corresponding to the larger shock velocity, while the 
cloud radius is in the range of (0.7-1.2) x 10 17 cm -2 with the 
lower limit corresponding to the lower shock velocity. If the 
cloud shape is clumpy, we come to the estimate of the cloud 
mass in the range of 10 -3 — 10~ 2 Af©. 

We also have performed simulations for the cloud B 
and find that the Ca 11 column density and acceleration are 
reproduced in cases of high and low cloud shock velocity for 
the cloud column density in the range of (4—5.4) x 10 17 cm -2 
and cloud number density in the range of 3.8-5.5 cm -3 , both 
values being similar to those for the cloud A. The fact that 
both clouds lying at the opposite sides of Vela SNR have 



comparable parameters indicates that low column density 
clouds may be rather common at least for the vicinity of 
the Vela SNR. We emphasise that this conclusion is rather 
speculative and requires confirmation by means of similar 
study of other high-velocity components in Vela SNR. 



5 DISCUSSION AND CONCLUSIONS 

Our goal was to study the interstellar absorptions in spectra 
of stars towards RX J0852. 0-4622 focusing primarily on the 
high velocity components associated with shocked clouds. 
We find that the low velocity absorption components do not 
show significant anomalies when compared to properties of 
absorbers in other surveys. The major differences with other 
surveys are probably caused by the selection effects due to 
different spectral resolution and S/N ratios. 

We find an interesting effect in the distribution of K I 
absorbers along the distance of background stars towards 
RX J0852. 0-4622: the lack of noticeable K 1 absorbers at the 
distances d < 600 pc followed by the "jump" of the strength 
of Kl absorbers for d > 600 pc which is indicative of the 
underpopulation of dense clouds at d < 600 pc. Generally, 
the distribution of H I gas in external spiral galaxies is inter- 
mittent with characteristic scale of several hundr ed parsecs, 
which is seen clearly, e.g., in M 51 21 cm map (|Rots et al.l 
Il990l ). It might well be that we look at the RX J0852. 0-4622 
through such local void with small number of cool clouds. 
Another possibility is that the hollow in the distribution of 
Kl absorbing clouds might be outcome of the Gum nebula 
dynamical evolution. Indeed, the Gum supershell is shown 
to be depopulated by ne utral gas clouds between 350-570 pc 
i|Woermann et al.ll200lT ). 

We have detected four high-velocity ( — v\ ST — > 
100 km s _1 ) Call clouds in spectra of three stars. Towards 
HD 75309 we see t wo shocked clouds, bot h observed earlier 
in 1993 and 1996 (|Cha fc Sembachll2000l ). The large time 
span of 12-15 yr permitted us to detect and measure the 
acceleration for cloud A and B which turn out to be com- 
parable (~ 10 -3 cm s -2 ). For the first time the acceleration 
of shocked clouds is used in combination with the velocity 
and Can column density to recover the cloud total hydro- 
gen column density (~ 6 x 10 17 cm -2 ), which implies the 
cloud radius of ~ 10 17 cm for the estimated cloud density 
n c o ~ 5 cm -3 . Regardless of our plane shock model is rather 
crude for the description of a blast wave/cloud interaction 
we believe that our model grasps the major features of this 
phenomenon. 

The fact that approaching and receding high-velocity 
clouds on the opposite sides of Vela SNR have comparable 
parameters indicates that they are possibly rather typical 
small scale (a ~ (1 — 2) X 10 17 cm) clouds in the vicinity 
of Vela SNR. It is not clear whether all the high-velocity 
clouds observed in Vela SNR have similar column density 
but, if we admit they do, it would be instructive to check 
the outcome of this conjecture. In our sample of 14 stars 
eleven do not show fast components which translates into 
the occultation optical depth (average number of clouds 
on the sight line) r oc = 0.24. This value is sampled at 
front and rear sides of the blast wave with the thickness 
A R « 0.03.R. Adopting the cloud number density distribu- 
tion dv Ida oc a~ 3 (cf. Appendix O which follows from the 
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observed proba bility density function of H I column density 
p(N) oc JV -1 (|Stanimirovic fe Heilesl [20051 ) and using the 
estimate r oc = 0.24 one finds the volume filling factor of 
clouds with the radii f < an < 2 to be / « 0.0f6. For 
the cloud total hydrogen density ~ 5 cm -3 the filling factor 
of the small size clouds extrapolated to larger sizes should 
not exceed 0.1 in order to be consistent with the average IS 
density of ~ 0.5 cm" 3 . This means that the upper limit of 
cloud size for that particular variety is ~ 4 x 10 17 cm. We 
therefore should admit that either the accelerated clouds are 
very rare and we catch an improbable chance, or, alterna- 
tively, the small-size clouds are rather numerous. If this is 
the case then possibly this variety of small clouds have not 
yet been observed. Indeed, if calcium were depleted in these 
clouds with 5 — 2 x 10~ 4 then the Can column density 
would not exceed 2 x 10 8 cm -2 which is below the detection 
limit. On the other hand, if Ca depletion we re moderate, 
e.g., 5 ~ 2 x 10" 2 likewise in clouds studied by I Bert in et al.l 
|l993t) then the Can column density would be as high as 
2 x 10 10 cm" 2 and these clouds could fall in the range of 
observed thinnest Can absorbers. 

At a first glance these small clouds could be identified 
with thin cl ouds (N s ~ 1.3 x 10 18 cm" 2 ) revealed by 21 cm 
absor ption (|Braun fe Ka nckar 2005]; IStanimirovic fc Heiled 
2005). The problem, however, is that these HI clouds are 
cool (T ~ 50 K) and therefore dense (n(H) ~ 30 cm" 3 ), 
factor ~ 5 more dense than our clouds. Another problem 
their low filling factor. Indeed, on the sight line towards 
3C286 extragalactic source (6 = 81°) only three clouds are 
detected, i.e., the average number of clouds on the sight line 
t oc =3. This value is translated into the filling factor / of 
clouds using the relation 

— !(*)'■ < 4 > 

where h z is HI scale height perpendicular to the Galac- 
tic plane, a is the cloud radius. For the typical density of 
~30 cm" 3 the cloud radius is a ~ 4 x 10 16 cm. Inserting 
into the equation Q the values of a, r oc , and h z = 10Q pc we 
get / « 5 x 10~ 4 , which is « 30 times smaller than the filling 
factor of the small clouds around Vela SNR derived assuming 
that all the high-velocity clouds are small-size clouds. This 
disparity and larger density of thin 21 cm absorbers imply, 
that the two varieties of small clouds seem to be different. 

Ultraviolet HST spectroscopy reveals in the direction 
(I, b) ~ (83° , —50°) a cloud with the column density o f JVh ~ 
6 x 1 17 cm" 2 and number density of 20-45 cm" 3 l|Welty| 
l2007t ). Again, this cloud has comparable column density with 
our A and B clouds but the number density is of one order 
magnitude larger. In fact, the density of A and B clouds (~ 
5 cm" 3 ) is somewhat unusual: it falls into unstable region of 
phase diagram for t he typical pressure P/k ~ 3000 cm" 3 K 
(|Womre et al.ll2003l ). This means that either A and B cloud 
were in the unstable state, or the pressure in the vicinity 
of Vela SNR is significantly (one order of magnitude) larger 
than the typical pressure. The latter, however, seems less 
likely because according to the extensive survey of UV fine- 
structure C I absorption lines the typical pressu re lies in the 
range 10 3 < P/k < 10 4 (| Jenkins fe TripplhoOll ). 

To conclude, the detection of the acceleration of two 
shocked clouds in Vela SNR along the same sight line and 
the fact of similarity of their parameters are exciting re- 



sults which might have interesting implications for our un- 
derstanding of cloudy structure of the ISM at small scales 
of ~ 10 17 cm. However, the issue, whether these clouds are 
common for diffuse ISM, or we observe highly improbable 
event, remains to be explored, e.g., by means of the study the 
velocity evolution of other high- velocity absorbing clouds in 
Vela SNR. 
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APPENDIX A: RELATION BETWEEN 
OBSERVED AND TRUE DISTRIBUTION OF 
CLOUD SIZES 

We consider homogeneous spatial distribution of spherical 
clouds of the same density. A line segment of a length L 
produced by a random sight line can be considered as a 
proxy for the corresponding column density N. The prob- 
ability density function of segments for a cloud of the ra- 
dius a is p(L,a) = (l/2)L/a 2 . Let the true cloud number 
density distribution of cloud radii be dv /da — Ca' 1 , while 
dn/dN oc dn/dL is the observed column density (or seg- 
ment) distribution per unit length along some line of sight. 
We then have an obvious relation between dv /da and dn/dL 

dn/dL — / 7ra 2 (dv / da)p(L, a)da oc L~ iri ~ ¥) , 

J L/2 

where we assume that 7 > 2 and a max 3> L. For the 
case dn/dN oc JV" 1 sugg ested by 21 cm observations 
l|Stanimirovic fc Hcilcs 20Q3) one gets 7 = 3. 
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